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Summary 

Neutrinos interact with matter only through weak processes with low cross-section. 
To detect cosmic neutrinos most efforts have relied on the detection of visible Vavilov- 
Cerenkov light in detectors embedded in the target volumes. To access the decreas- 
ing flux of ultra-high energy neutrinos, far above 1 PeV, ideas on how to increase 
the detection volume by observing coherent radio frequency emission caused by the 
Askaryan effect have been put forward. Here we describe how a satellite in lunar 
orbit equipped with an electromagnetic vector sensor could detect Askaryan pulses 
induced by neutrinos interacting with the moon. The threshold neutrino energy is 
found to be 5 x 10 19 eV for this setup, and the sensitivity is determined from simu- 
lations. A model dependent event rate of 2.2 events per year is calculated. 



Introduction 



The high end of the cosmic neutrino energy spectrum remains unstaked territory. 
Different models (Active Galactic Nuclei 01 , Topological Defects 01, GZK-cutoff 
JH@] ) predict sources of cosmic neutrinos with energy far above 1 PeV (10 15 eV), 
termed ultra-high energy (UHE) neutrinos. Validation of these models remains an 
experimental challenge. 

Detection of cosmic neutrinos is possible by observing optical Vavilov-Cerenkov 
emission from secondary particle showers produced in dense media from neutrino 
interactions. With the expected decrease of the neutrino flux with energy |5], larger 
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detector volumes are necessary for higher energy particles. The ICECUBE detector 
will equip 1 km 3 |6J] of South pole ice with photomultiplier tubes and still larger 
volumes will be needed to access the UHE neutrinos. It becomes increasingly im- 
practical to detect the UHE neutrinos by their induced optical radiation, also since 
the optical output is generated in an incoherent process and therefore the output 
power increases only linearly with the neutrino energy. 

One way to overcome this limitation and access a sufficient detector volume is to 
detect Vavilov-Cerenkov radio emissions. At radio frequencies, the emission is gen- 
erated coherently and the output power is proportional to the square of the primary 
energy. Coherent emission is possible since the showers develop a negative charge 
excess going through the material; this is known as the Askaryan effect 0], j^]. 

Radio methods for neutrino detection have been employed by several groups world- 
wide. One such experiment was performed by the GLUE collaboration, who utilised 
the Goldstone radio telescopes in search for pulses from the Moon |9]. Embedded 
aerials in the Antarctic ice have been used by the RICE 1 10] group. Data from the 
FORTE satellite II ill has been searched for pulses generated in the Greenlandic ice 
sheet. No events were reported from any of these experiments, and this sets upper 
limits on the neutrino flux above 10 17 eV. Of great promise for the future is ANITA 
fi3l> a balloon-carried transient detector to be launched in 2006. 

In this report, we will discuss how a digital broadband vector radio receiver on board 
a satellite orbiting the Moon could be used to detect radio pulses from neutrinos 
interacting with the lunar surface. To detect Askaryan pulses from the Moon by 
satellite measurements has previously been suggested by Gorham fTlll . 

Radio emission from particle showers 

Neutrinos interact with matter only through weak processes with low cross-sections. 
The cross-section increases slowly with energy, and determines a minimum neutrino 
energy required for interaction with the lunar material. This energy can be estimated 
by comparing the energy dependent interaction length, given by the cross-section, 
to the lunar radius = 1738 km. Using an average density of 3.34 g/cm 3 for the 
lunar interior 1 14J gives a value of E m { a ~ 3 x 10 16 eV. Above this energy, the Moon 
is opaque to neutrinos. 

After the first interaction, a shower of secondary elementary particles is gener- 
ated. The charged component consists mainly of a compact region of electrons and 
positrons propagating at superluminal speeds 1 15]. Perpendicular to the direction of 
motion, the spatial scale of the shower is given by the Moliere radius Q(J, which 
for Si is about 5 cm ll5J . On the average, 90% of the shower energy is contained 
inside a cylinder of this radius around the shower core. 
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Radiation in homogeneous media 

For a point charge moving with uniform velocity in a homogeneous, dielectric 
medium of refractive index n = -y/e, we can obtain, in the standard way, the Fourier 
transform of the electric field by direct solution of the Maxwell's equations fl7ll . 
Considering satellite detection, Fraunhofer geometry applies. The vectorial dis- 
tance traversed by the charge is denoted d, with a component dx perpendicular to 
the direction to the observation point. Following 111 811 . the radiated electric field 
becomes 



sin0 







(1) 



Eo , (x) = ^W^ e i« B -M. dxe 
4neoc 2 r 

where the square-bracketed expression resembles the angular dependence of a slit- 
diffraction pattern. An overall phase is provided by t\, the time when the particle 
starts to radiate. The phase-factor 

ltd 

<j> = — (1-njBcosfl), (2) 

is defined in terms of d — |d|, the observation angle from the shower axis, and the 
velocity factor /3 = \\\/c. If the particle velocity exceeds the phase speed of light 
in the medium v = c/n, there is a maximum field strength attained for § = 0. This 
condition defines the Cerenkov angle 6 C = arccos l/(«j3). When the phase factor is 
small, a McLaurin expansion of the phase dependence in is possible. Retaining 
only the first term, the field expression simplifies to 

E4x)-^^- k ^d ± . (3) 

The necessary criterion for (|3j to be valid is either that the wavelength is sufficiently 
long (d <C X) or that the observation angle is close to the Cerenkov angle. 

The total field from a shower is obtained by adding the contributions of all particles. 
The result depends on the spatial scales, as compared to the wavelength of the ra- 
diation. In the optical regime, all particles radiate incoherently and the individual 
phases must be taken into account through Q. The individual phase factors are uni- 
formly distributed which makes the signs of the charges unimportant and the field 
strength scales with the square root of the primary energy. 

The radio regime is at the other end of the spectrum. If all spatial scales of the 
shower are assumed to be smaller than the wavelength, the radiation will be emit- 
ted almost coherently and the field strength will scale linearly with energy. In this 
regime, Eq. (|3j is valid, and only the net charge gives rise to the field. 

Between the two regimes there is a transition regime of partial coherence. Closer 
to the Cerenkov angle, higher frequency components display coherence. Following 
llisil . we define the highest frequency with coherent emission at the Cerenkov angle 
as the specific decoherence frequency of the material. 
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The development of charge asymmetry 

For a shower with equal numbers of electrons and positrons, no coherent radiation 
can be emitted. Askaryan realized that when such a neutral shower develops in 
a material, scattering processes will create an imbalance with a negative charge 
excess 0], The most important contribution comes from Compton scattering (7 + 
e atom ~~ * 7+ e ~) which incorporates new electrons into the shower. Another process 
increasing the number of shower electrons is Bhaba scattering (e + + e~ tom — > e + + 
e~). Finally, there is annihilation (e + +e~ tom — ► 7+ 7), further reducing the positron 
component. Typical values for the fully developed charge excess are about 20% 
Q. 

Spectral parametrization 

The total radio pulse spectrum has been determined from simulations taking into 
account the primary neutrino interaction and the following shower development. 
Such simulations for showers in ice llisll of total energy Wo have been generalized 
to other dense media and can be summarized in the empirical formula 11911 

i E ^<«' = 7 (nsv) f 1+<M ^)-" [Vn, " MH "' 1 ■ (4) 

The amplitude A is determined empirically, as is the specific decoherence frequency 
COn. The general features of this formula have been verified experimentally at SLAC 
12(11 . With bunches of GeV photons released in a Si target, a total shower energy 
of 10 19 eV was attained. We will adopt the experimental parameters for Si to the 
radiation from the Moon, namely A = 2.53 • 1CT 7 V MHz -1 and C0q = 2.5 GHz. 

Showers from primaries at an energy above 1 PeV interacting initially through the 
electromagnetic channel are known to be elongated by the Landau-Pomeranchuk- 
Migdal (LPM) effect I21I l22ll. This significantly reduces the angular spread of 
the Vavilov-Cerenkov cone 11911 . making such showers unfavourable for detection. 
They will therefore not be considered in the following. Neutrinos interacting ini- 
tially through the hadronic channel will not suffer LPM elongation since many par- 
ticles of lower energy are generated at an early stage of the shower development. 
An average fraction 0.2 of the primary energy is deposited in the hadronic shower 
component 1 23] . 

The coherent radiation from showers of many particles is not limited to the Cerenkov 
angle, but, as discussed above, for low frequencies the coherence is maintained 
further off-axis. The reduction in field strength off-centre approximately follows a 
Gaussian: 

|E fi) (r,0)| = |E (B (r,0 c )| e -K^) 2 ( 5 ) 
with a frequency dependent width given by AO w 2.4° C0o/ CO. 
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The Moon as a neutrino detector 



The uppermost layer of the lunar surface consists of a fine-grained material known 
as the regolith. The electrically neutral regolith is composed of single mineral 
grains, small rock fragments and combinations of these, fused by impacts. The 
depth of the regolith is believed to be at least 2-8 metres in the maria, increasing to 
15 meters or more in highland regions 1 14]. Silicates are the most abundant com- 
ponent in the regolith, but important for the dielectric characteristics are also the 
varying concentrations of Fe and Ti oxides. 

Assuming the mass density of the regolith layer to be constant p = 1 .7 g/cm 3 over 
the depths considered (0-10 m), the dielectric constant is taken to have an average 
value of e = n 2 = 3. It has been found previously that the value of the loss tangent 
depends most significantly on the concentrations of TiC>2 and FeO in the material. 
Assuming an average total concentration of these metallic oxides of 5%, the loss 
tangent tan5 = 0.003 1 24] . This gives an attenuation length of 100 m for a 100 
MHz signal. The attenuation length is inversely proportional to the frequency. 



Satellite detection of Askaryan radio pulses 



Threshold 



From the empirical formula © it is possible to estimate the threshold energy for de- 
tection, assuming the background level is known. The required signal level is com- 
monly given as a power figure in dBm per Hz, related to the spectral field strength 
through 

P(dBm) = 101og 10 (T^pr) , (6) 

1 q6 I PR 

E <» = ^\^ [juVm^MHz" 1 ] (7) 
flA V Aa> 

where /ia is the effective antenna height, R = 50 D. and A© the system bandwidth. 
Equating with © now gives the threshold energy 

This analysis is valid in the regime where E m increases linearly with co [cf. Eq (4)], 
i.e. for (O + Aco <C coo- 



Aperture and event rate 



In order to estimate the expected event rate, a calculation of the effective lunar aper- 
ture (area x solid angle) for neutrino detection as a function of energy is required. 
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Figure 1: Effective aperture for lunar neutrino detection. One sigma error bars 
indicate statistical deviations. 



The expected event rate is then obtained by convolving the differential flux with the 
aperture 

Where the aperture AQ. is factorized as the geometric aperture [Ai2] geom multiplying 
a specific efficiency, S. We determine this efficiency by Monte Carlo simulations. 

Since the angle for total internal reflection at the Moon-vacuum interface 0xir is 
complementary to C , only emission from showers that are upgoing with respect to 
the surface will escape when the Cerenkov cone is narrow. For lower frequencies 
the cone is wider, and total internal reflection is not a significant problem since the 
transmission rises to a maximum just a few degrees off &tir- 

All rays will be refracted at the interface. For a planar surface, the refraction will 
follow Snell's law. This will broaden the Cerenkov peak and the acceptance angle 
is increased on the expense of the field strength. A more elaborate simulation of the 
interface transmission should take the structure of the surface into account. 

Experimental equipment 

The ELVIS (ELectromagnetic Vector Information Sensor) instrument l25ll recently 
proposed for the Indian Moon mission Chandrayaan-1, scheduled for launch in 
2007-2008, provides a suitable test-bed for detection of lunar neutrinos by radio 
methods. The equipment consists of three 5 m dipole aerials in an orthogonal con- 
figuration, providing full 3D sampling of the electric field vector and the advantage 
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Figure 2: Aperture in the interesting energy regime close to threshold with a poly- 
nomial fit to the simulation results. 



of polarimetric measurements. The signals are digitized in linear 14 bit ADCs. The 
sampling frequency is currently limited to 80-125 Msamples/s, but when operated 
in transient mode, signals up to 200-500 MHz are detectable. 

Simulation results 

Calculations of the aperture and expected event rate has been performed for the 
ELVIS and Chandrayaan- 1 setup using our own Monte Carlo code. The detector 
parameters used were 100 km height above the surface, CO — 100 MHz for the obser- 
vation frequency with an estimated bandwidth of 50 MHz. The minimum detectable 
signal strength was taken as a conservative —135 dBm. The threshold according to 
formula (|8j becomes 5 x 10 19 eV, but, as can be seen from the results, the aperture 
below 10 20 eV is effectively negligible. 

The aperture was calculated using cross-sections from Gandhi et al. |26j, ZHS 
parametrization of the field according to Q, full ray-tracing and smooth surface in- 
teractions as outlined above. The result is presented in Figure 1, showing a compli- 
cated energy dependence near threshold followed by a logarithmic increase of aper- 
ture with primary energy. The region close to threshold is presented more clearly in 
Figure 2, with a polynomial fit used for event rate integration. 

Adapting the GZK neutrino flux of Engel et al. i27Tl as a standard neutrino source at 
ultra-high energy, the integrated event rate over a year of observation time becomes 
2.2 expected events for these parameter values. 
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Conclusions and outlook 

Radio detection of coherent Vavilov-Cerenkov emission originating from showers 
induced in the moon by ultra-high energy cosmic neutrinos has been reviewed. Cal- 
culations for a satellite detection experiment reveals the very low event rate to ex- 
pect, but it should be possible to improve the event rate by orders of magnitude 
through optimization of the experimental parameters. 

Highly interesting is also to study the use of an antenna array on the lunar surface 
for neutrino detection. This must be performed through extensive simulations. Such 
simulations are underway, giving further results and more insights in the near future. 
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